Is there hydrodynamic flow at RHIC ? 
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It is argued that the observation of anisotropic azimuthal distribution of final state particles 
alone is insufficient to show whether the formed matter at RHIC behaves like hydrodynamic flow. 
Examining the intrinsic interaction (or correlation) of the formed matter should provide more definite 
judgement. To the end, a spatial-dependent azimuthal multiplicity-correlation pattern is suggested. 
It shows clearly in the pattern that there are two kinds of interactions at the early stage of Au -|- 
Au collisions at ^snn ~ 200 GeV generated by RQMD with hadron re-scattering and AMPT with 
string melting. This is out of the expectation from the elliptic flow driven by anisotropic expansion. 
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The data from current relativistic heavy ion experi- 
ments show that a new form of matter — quark-gluon 
plasma (QGP) has been produced at RHIC [H, 0. The 
second Fourier coefficient V2 of the anisotropic transverse- 
momentum pt distribution of final state particles is be- 
lieved to provide the anisotropic collective flow behavior 
at the early stage of collision. The successful hydrody- 
namic description [3] on the observed mass dependence 
of V2 at pt < 2 GeV shows that the observed dense mat- 
ter behaves like a perfect fluid rather than an ideal gas, 
and is, therefore, referred to as sQGP. 

However, hydrodynamics can still not quantitatively fit 
the observed mass dependence of V2 P, Q ■ The recently 
measured elliptic flow V2 from Cu -|-Cu collisions at 200 
GeV is unexpected as large as that from Au -I- Au colli- 
sions at the same energy [^,@. Moreover, the resulting 
matter may be treated as hydrodynamic flow only if the 
initial interaction among the constituents are sufficiently 
strong to establish local thermal equilibrium rapidly, and 
then to maintain it over a significant evolution time. No 
known strong-interaction process could be thermalized 
on such a short timescale. A Liquid without viscosity is 
also hard to be understood theoretically 0]. So, there ap- 
pear a number of alternative non-equilibrium treatments, 
which have also been compared to RHIC data [1, ■ 

To conclusively clarify the debate, a direct experimen- 
tal examination on the intrinsic interaction of the formed 
matter is neccessary. The hydrodynamic flow at RHIC is 
supposed to be driven by the so called anisotropic expan- 
sion. In non-central collisions, the initial participant zone 
of the two colliding nuclei is approximately an ellipse, 
and the density gradient along the short side of ellipse is 
larger than that along the long side. It is argued that the 
larger density, or pressure, gradient along the short side 
of ellipse makes collective expansion to be privileged in 
this direction, i.e., the anisotropic expansion, producing 
in-plane elliptic flow, or the transverse-momentum of fi- 
nal state particles distribute in an ellipse perpendicular 
to the one in coordinate space. 

However, the main physical quantity which can be ex- 
tracted from experimental data in exploiting relativis- 
tic hydrodynamic approach is the elliptic flow parameter 
V2 ■ It only indicates the possible preferential direction of 



expansion and contains no information on the intrinsic 
interaction of the formed matter. It therefore is insuffi- 
cient to assure whether the formed matter behaves like 
hydrodynamic flow. 

If the anisotropic expansion is the only driver of the 
elliptic flow, the distribution of intrinsic interaction (or 
correlation) of flow should have the same anisotropy, i.e., 
in-plane like. Moreover, if it is really hydrodynamic flow, 
it should be well locally thermalized and reach thermal 
equilibrium. Then all other interaction history before 
anisotropic expansion should be forgotten. These char- 
acteristics can be examined in an experimentally measur- 
able correlation pattern. 

In this letter, we will first introduce the spatially- 
dependent correlation pattern, i.e., neighboring angular- 
bin multiplicity correlation pattern. Then, we demon- 
strate that at least two kind of interactions are re- 
vealed by the suggested correlation pattern in Au -I- 
Au colhsions at 200GeV, generated by RQMD [3 and 
AMPT [l^. Finally, how to experimentally measure the 
correlation pattern and anisotropic correlation coefficient 
is discussed. 

To examine the intrinsic interaction of highly 
anisotropic system, a spatial-dependent bin-bin correla- 
tion is called for. Conventionally, the spatially averaged 
bin-bin correlation has been used in multiparticle pro- 
duction in exploring self-similar fractality , where the 
system is supposed to be homogeneous, and only scaling 
in the shrinking of phase space is concerned. Another 
intrinsic interaction related measure is the 2-particle az- 
imuthal correlation 0, [l^l ■ It concerns the average cor- 
relation of two particles separated by a certain angle, 
no matter where the two particles are in the azimuthal 
space. It therefore can not tell us where the preferential 
direction of intrinsic interactions are. 

The newly suggested spatial- dependent neighboring bin 
correlation pattern [l^ provides a typical spatial distri- 
bution of two-bin correlation. The information on intrin- 
sic correlation can be well presented by the measure, and 
it should give more direct and definite judgement on the 
properpty of the formed matter at RHIC. 

It is well-known that the general 2-bin correlation is 
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defined as 
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where mi and m2 are the positions of the two bins in 
phase space and rim is the measured content in the mth 
bin. 

We divide the 2tt azimuthal angle equally into M bins 
and specify the multiplicity in the mth angular bin. 



If we let mi = m and m2 = m + 1, C„ 



is reduced to 



the neighboring angular-bin multiplicity correlation pat- 
tern. 
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It is clear that the correlation pattern measures how 
the nearby particles correlate with each other in differ- 
ent directions of azimuthal space. If the particles arc 
produced independently in the whole phase space, then 
{nynnm+i) = {nm){nyn+i), and Crm^ra^ vanishes. 

In order to apply this correlation pattern to current 
relativistic heavy ion collision, we choose the RQMD 
and AMPT models as examples. The RQMD (relativis- 
tic quantum molecular dynamics) with rc-scattcring is 
a hadron-based transport model [15| . The final hadron 
interactions are implemented in the model by hadron re- 
scattering. The anisotropic collective flow produced by 
the model is much smaller than the observed data at 
RHIC. In contrary to the RQMD model, the AMPT is 
a multi-phase transport model, where both hadron and 
parton interactions arc taken into account. In the AMPT 
with string melting, the parton level transport is fully 
taken into account, and the observed anisotropic collec- 
tive flow at RHIC is well reproduced 

For Au -I- Au collisions at -^/snn = 200 GeV, we gener- 
ate 249,824 and 204,004 events using RQMD with hadron 
re-scattering and AMPT with string melting, respec- 
tively. Their neighboring angular-bin multiplicity cor- 
relation patterns are shown in Fig. 1(a) by open and 
solid circles, respectively. Here we partition the whole 
azimuthal range 27r uniformly into 50 equal size angu- 
lar bins. = refers to the direction of the reaction 
plane in nuclear collision. The errors are statistical only 
and most of them arc smaller than the symbol size in this 
and following figures. It is clearly shown in Fig. 1(a) that 
correlation patterns from these two models are — cos 2(f) 
(out-of-plane) like, opposite to the well-known cos 2(j) (in- 
plane) liked azimuthal distribution. This is in contrary 
to the expectation that the formed matter expands col- 
lectively toward in-plane direction. Some unexpected in- 
teractions should be responsible for such a result. 

In order to see how the results come, the centrality 
dependence of neighboring angular-bin multiplicity cor- 
relation patterns from these two models are presented in 
Fig. 1(b) and (c), respectively, where three typical cen- 
tralities are specified in the legends. One can observe that 
the two models give qualitatively the same centrality de- 
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FIG. 1: (a)The neighboring angular-bin correlation patterns 
for Au 4- Au collisions at 200 GeV from the RQMD with 
re-scattering and AMPT with string melting. The centrality 
dependence of the correlation patterns from (b)the RQMD 
with re-scattering , and (c) AMPT with string melting. 



pendence of azimuthal correlation pattern. The correla- 
tion patterns are cos 20 like in peripheral collisions, then 
turn to flat in mid-central collisions, and become — cos 20 
like in near-central collisions. Here, we present only three 
centrality ranges to show their typical behavior. In fact, 
the correlation pattern changes gradually from cos 20 to 
— cos 20 with centrality. It is clear that two opposite 
trends dominate in peripheral and near-central collisions, 
respectively. In the mid-central collisions, the two trends 
turn to balance and the correlations become equal in all 
directions. Moreover, these characteristics are indepen- 
dent of the specific assumptions implemented in the two 
models, in particular independent of the hadronization 
schemes assumed in the models. 

The characteristics of correlation pattern reveal that 
there are two opposite intrinsic interactions in the formed 
matter in these two transport models. One has the same 
preferential direction as the anisotropic expansion. The 
other one is opposite to it. This also shows that the 
anisotropic azimuthal distribution is not only driven by 
anisotropic expansion. It is resulted from the combina- 
tion of these two opposite interactions. 

The anisotropic expansion and the late hadronization 
are impossible to produce strong correlations in out-of- 
plane direction. Only the initial source eccentricity in 
non-central collisions is preferential in the direction. It 
results in a larger initial number of participant nucleons 
in the out-of-plane direction, which in turn could gener- 
ate stronger interaction in the direction. As long as the 
system is not fully thermalized, this initial interaction 
will compete with the subsequent anisotropic expansion. 

In peripheral collisions, the overlap zone is small and 
so is the number of participant nucleons, but the differ- 
ence between the minor and major axes of overlap ellipse 
is large, and so is the difference of pressure gradients. In 
this case the anisotropic expansion dominates the final 
observables, and the effects of initial interaction in corre- 
lation patterns are hidden. In near-central collisions, the 
overlap zone becomes large and the difference between 
minor and major axes of ellipse is small, so that the ini- 
tial interactions are strong enough to show themselves up 
in final observable. This is why the out-of-plane correla- 
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FIG. 2: (a) Transverse-momentum and (b) rapidity depen- 
dence of U2 at three centralities for Au -I- Au collision at 
200GeV from RQMD with re-scattering. 



tion patterns appear at near-central collisions. 

So the behavior of the formed matter in these two 
transport models are far from the flow in relativistic hy- 
drodynamics sense. This is out of the current expectation 
for the formed matter at RHIC. Measuring the correla- 
tion pattern by the data of relativistic heavy ion collisions 
at RHIC and LHC is therefore looking forward. As long 
as the observed preferential direction of correlation pat- 
tern arc different from that of its azimuthal distribution, 
such as what we show by transport models, then there 
should be no hydrodynamic flow at RHIC. On the other 
hand, if the experimentally measured correlation pattern 
has the same anisotropy as its azimuthal distribution, it 
will be a strong support to the current expectation at in- 
trinsic interaction level that the formed matter at RHIC 
indeed behaves like hydrodynamic flow. 

The correlation patterns are typical periodic functions 
of azimuthal angle in peripheral and central collisions, as 
shown in Fig. 1. So they can be well expanded by Fourier 
series. 



0,0+50 



c 



1 + 2Ui COs(l(0 — Ipr)) 



(3) 



where the V'r is the direction of reaction plane, and is 
zero in the model analysis. But in real experimental 
data analysis, it has to be determined event-by-event, 
and thereby refers to event-plane. It has been carefully 
estimated in the measurement of anisotropic elliptic flow 
V2 in current relativistic heavy experiments p^ . [l9| . The 
main contribution in the expansion series comes from 
cos 2(0 — V-'r)- Its coefficient U2 provides the preferen- 
tial direction and strength of anisotropic correlation pat- 
tern. We specify it as anisotropic correlation coefficient 
(ACC). It will make the systematic study of the correla- 
tion pattern easy. 

It is interesting to see how ACC, U2, depends on the 
transverse-momentum pt of final state particles, in com- 
parison to the corresponding pt dependence of V2. It 
is known that the evolution schemes of RQMD with re- 
scattering and AMPT with string are different. In the 
former, the pt spectrum is determined by the temperature 



of thermal source. High pt particles are emitted early at 
high temperature and low pt ones are emitted later on at 
low temperature. The range of pt of final state particles 
is related to its emitting proper-time [l^. So the pt de- 
pendence of U2 in RQMD with re-scattering will present 
how the correlation pattern changes with evolution. 

The results are presented in Fig. 2(a). We can see that 
in each pt interval, U2 keeps positive in peripheral col- 
lisions, becomes negative for mid-central collisions, and 
becomes even more negative for central collisions. They 
are similar to that for all pt particles shown in Fig. 1(b). 
It should also be noticed in Fig. 2(a) that U2 is almost 
independent of the choice of pt ranges of final state parti- 
cles in peripheral and mid-central collisions, but decrease 
rapidly with the increase of pt in central collision. This is 
understandable since highpt particles arc emitted earlier, 
and less influenced by the later anisotropic expansion in 
central collisions. 

On the contrary, each parton in the AMPT with string 
melting has its own freeze-out time, which span a long 
period after the initial interaction of the two nuclei, and 
are unrelated to each parton's transverse momentum [igI . 
[iTj . So similar pt dependence of U2 can be observed in 
this model only when the chosen interval of pt is very 
large. 

The rapidity dependence of ACC, U2, is further stud- 
ied by these two transport models. They give qulitativcly 
the same dependency. The results from RQMD are pre- 
sented in Fig. 2(b), where three typical rapidity ranges, 
i.e., forward, backward and central rapidity ranges, are 
chosen. It shows that the correlation pattern is indepen- 
dent of the choice of rapidity range, and similar to that 
in the whole rapidity space. So in finite rapidity ranges 
of current relativistic heavy ion experiments [l9| , study- 
ing the correlation pattern and anisotropic correlation 
coefficient is expectable, and will provide more definite 
evidence on whether or not the formed matter behaves 
like hydrodynamic flow. 

To the summary, it is argued that the observation of 
anisotropic azimuthal distribution of final state particles 
alone is insufficient to assure whether the formed matter 
at RHIC behaves like hydrodynamic flow. Examining the 
intrinsic interaction (or correlation) of the formed mat- 
ter should provide more definite judgement. To the end, 
a spatially-dependent azimuthal multiplicity-correlation 
pattern is suggested. It shows clearly that there are two 
kinds of interactions at early stage of Au -I- Au collisions 
at y'sNN = 200 GeV, generated by RQMD with hadron 
re-scattering and AMPT with string melting. One is in- 
plane preferential as expected from anisotropic expan- 
sion due to initial eccentricity in non-central collisions. 
Another new one is out-of-plane preferential, which may 
be resulted from the larger initial number of participant 
nucleons in these direction. These characters of correla- 
tion pattern show at least in two transport models that 
the formed matter does not behave like hydrodynamic 
flow, in contrary to current expectation. Finally, how 
to experimentally measure the correlation pattern and 
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anisotropic correlation coefficient is discussed. 
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